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We study an effect of inhomogeneity of density distribution of the Universe. We propose 
a new Lagrangian perturbation theory with a backreaction effect by inhomogeneity. The 
inhomogeneity affects the expansion rate in a local domain and its own growing rate. We 
numerically analyze a one-dimensional plane-symmetric model, and calculate the probabil- 
ity distribution functions (PDFs) of several observed variables to discuss those statistical 
properties. We find that the PDF of pairwise peculiar velocity shows an effective difference 
from the conventional Lagrangian approach, i.e. even in one- dimensional plane symmetric 
case, the PDF approaches an exponential form in a small relative-velocity region, which 
agree with the N-body simulation. 

I. INTRODUCTION 

The present Universe shows a variety of structures. How such a structure is formed in the evolution of the 
Universe. One of the most plausible explanations is that the nonlinear dynamics of a self-gravitating system 
provides such a scale-free structure during the evolution of the Universe. In order to clarify whether such a 
dynamics really gives an appropriate observed feature, we have so far three approaches: iV-body simulation, 
the Eulerian perturbation approach and the Lagrangian one. Although the final answer for a structure 
formation would be obtained by the iV-body simulation, it may be difficult to obtain enough resolution to 
discuss a fine structure of the Universe such as a scaling property. As for a perturbation approach, however, 
it is just an approximation and will break down in a nonlinear regime, although the Lagrangian approach 
would be better if we are interested in density perturbations. This is just because a density fluctuation 6 and 
a peculiar velocity v are perturbed quantities in the Eulerian approach , while a displacement of particles 
from uniform distribution is assumed to be small in the Lagrangian approach [^-^|. Its first order solution 
is the so-called Zel'dovich approximation (ZA) [0. The Lagrangian approach is confirmed to be better than 
the Eulerian approach by comparison of exact solutions in several cases . Therefore, we will adopt the 
Lagrangian approximation in this paper and discuss about how to improve it. 

In the standard approach of Newtonian cosmology, the global cosmological parameters such as Hubble 
expansion rate and mean density are given first by a solution of the Einstein equations, i.e. the so-called 
Friedmann-Robertson- Walker (FRW) universe, which is an isotropic and homogeneous spacetime. According 
to observation, however, a local structure in the Universe is definitely not homogeneous and isotropic pj. 
In the standard approach, the density averaged over the whole space (or a horizon scale) is assumed to be 
the energy density of the FRW spacetime. However, here the problems of how to average inhomogeneous 
matter fluid and how to define an averaged isotropic and homogeneous spacetime are arisen. This problem 
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has been discussed in the framework of Newtonian cosmology and general relativity by many authors |llj |l7j] . 
The most important point to discuss is how inhomogeneities affect the expansion law of the Universe and 
whether such a backreaction affects the evolution of the density perturbations. This is our subject to discuss 
in this paper. Although many authors proposed different methods of averaging procedure in general relativity 
JTT| p"5[ , a serious problem remains. Different gauge choices make the problem more complicated. 

Since a relativistic effect may not be so important, in this paper, in order to avoid such a difficulty, we 
discuss only an averaging procedure in the Newtonian framework. Proposing a averaging procedure, which 
is defined by spatial average of physical quantities, Buchert and Ehlers lead the averaged Raychaudhuri's 
equation [[L6| . The equation describes how the averaged expansion rate of domain with a finite volume 
evolves. This equation has an additional term, which we call a 'backreaction term' of inhomogeneities 
on averaged expansion. Then, Buchert, Kerscher, and Sicka estimated the backreaction term using the 
conventional Lagrangian perturbation approach ]lg| ] . They first consider density perturbations in Einstein- 
de Sitter (E-dS) Universe, then calculated the backreaction term and solved the averaged Raychaudhuri's 
equation. They showed difference between a cosmological parameter such as Hubble expansion rate in their 
averaged model and that of the E-dS model. Although they included a backreaction term to estimate the 
averaged variables in a local domain, they used the perturbed quantities from the E-dS Universe. In other 
words, they have not take into account a backreaction on the evolution of density perturbations. 

Here we improve their approach, i.e. wc include a backreaction effect to averaged expansion rate and 
solve the averaged Raychaudhuri's equation (the generalized Friedmann's equation) with evolution equation 
of perturbations in the averaged domain. In an averaged domain, the averaged density is either higher or 
lower than that of the E-dS universe. This difference will change the evolution of density perturbations. In 
fact, if the domain is overdense, growth rate of perturbations behaves as that in the closed universe, While, 
if it is underdense, it is just like a solution in the open universe. 

This paper is organized as follows: In §2, we shortly derive the generalized Friedmann equation, following 
Ehlers and Buchert, and estimate a backreaction term. We present our formalism in §3. We show some 
relation between two Lagrangian descriptions (the conventional one and ours) and give our initial setting in §. 
4. In §5, we analyze a simple example, i.e. a plane-symmetric one-dimensional model, to show the probability 
distribution of the observed quantities such as Hubble parameter, or density fluctuation, or peculiar velocity. 
The conclusion and discussion follow in §6. 



In the Newtonian cosmology, the expansion of a domain is influenced by inhomogeneity inside the domain. 
Such an effect may be evaluated by spatial integration of field variables in the Lagrangian domain, which 
evolves with matter fluid. Hence in this paper we study fields averaged over a simply-connected spatial 
Lagrangian domain T> at time i, which evolved out of the initial domain T>i at time i^. The locally averaged 
scale factor ap, depending on the content, shape and location of the domain T>, is defined by the volume of 
domain Vx>(t) = \T>\ and its initial volume = \T>i\ as 



We define a spatial averaging for any rank tensor field A = {Aij. (r, t)} by the volume integral normalized 
by the volume of the domain as 



II. AVERAGING OF INHOMOGENEITY 



A. The generalized Friedmann equation 




(2.1) 




(2.2) 



The average density is then given by 
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where the total mass Mjy = j v d 3 r p(r, t) for a domain T> is conserved. Using this averaging, we can derive 
the generalized Friedmann equation |lq| , which we will shortly derive here. 
We start from the hydrodynamic equations for a self-gravitating perfect fluid; 

^ + V r ■ (pu) = , (2.4) 

-^ + (u-V r )u = --V r P-V r $, (2.5) 

V r xg = 0, (2.6) 
V r • g = A - AnGp . (2.7) 

Although it is easy to extend the case with pressure of fluid, we consider just dust matter and ignore a 
pressure term. 

To discuss inhomogeneity, a spatial derivative of field will play a crucial role. In particular, it turns out 
that the spatial derivative of velocity u, which is divided into three variables as 

L . 

U hj = (7y + -OijV + OJij , 



is very important, where 



with Uij = -e ljk uj k , (2.8) 



6 = V • u (expansion) 
cry = + Uj ti ) - ~5ij6 (shear) 

u = -V x u (rotation). (2.9) 



The magnitudes of shear a and of rotation u> are defined by 



-(TijCTij ) , W = \u\ = ( -WijWij ) ■ (2.10) 



Using these variables, we find the basic equations for density, shear and rotation from Eqs. (|2.4|), (|2.5|) as 



p = -0p, (2.11) 
2 

us = --6u) + a-u> , (2.12) 

6 = A - AirGp - \d 2 + 2{uj 2 - a 2 ) . (2.13) 

The volume of a spatial domain T> described by its initial domain T>i through a transformation from the 
Eulerian coordinates r to the Lagrangian one q as 

V v = [ d 3 r= [ d 3 qJ r , (2.14) 

JV JVi 

where J r is the Jacobian of the transformation 

J r =det f|^J . (2.15) 

The change rate of the volume of a domain is then given as 

1 dVv 1 d f . T _ 1 f ~ 



d q J r = 77- J r 



Vd dt VvdtJ v . V v 



= — I d 3 q 6J r = — I d 3 r 9 = (6) v , (2.16) 
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where d/ 'At denotes the Lagrangian time derivative. Since {9) v is the averaged volume expansion rate of the 
domain T>, it may be natural to define the effective Hubble expansion rate of the domain T> by -ffp = (9)^ /3. 
From the definition of the volume Vp, we find Hx> = ap/a-p. 

The Lagrangian time derivative does not commute with spatial averaging. For an arbitrary tensor field 
A, Ehlcrs and Buchert showed the commutation rule: 



d (p)v 



(2.17) 



Applying this commutation rule to the basic equations, we find 

= -{9) v (p) v , (2.18) 



= - (9) v (u,) v + (« • V r u> D , (2.19) 
= A - AkG ( P ) v + | (0% - + 2 ((c 2 ^ - <a 2 ) p ) . (2.20) 



Replacing (#)•£> in Eq. (2.20) with 3du/ap , we find the averaged Raychaudhuri's equation 

3— + 47rG (p) v -A = Qv, (2.21) 



^4((a-(^)+ 2 (a-a) (2.22) 



where 

2 
3 

describes the backreaction due to inhomogeneity of the universe. For Qx> = case, the equation becomes 
the conventional Friedmann equation with a scale factor a(t) for homogeneous and isotropic universe. This 
is regarded as the generalized Friedmann equation. 



B. Evaluation of Backreaction Term Qx> 

In order to evaluate the back reaction term Qt>, it may be convenient to introduce three principal scalar 
invariants as follows: For 2-rank tensor field A — (Ay) in Cartesian coordinates, those are defined by 

I(A ij )=tv(A ij ) , (2.23) 

II(A y ) = i(tr(^) 2 -tr((A y -) 2 )) , (2.24) 

III(Ay) = det(A y ) . (2.25) 
In particular, for the velocity gradient of matter fluid (uj,j) , we find 

I(ttij) =Ui,i = V r .u = 0, (2.26) 

n K,i) = 2 {( u m) 2 ~ "••.-",••<} 

= iv r • {u(V r • u) - (u • V r )u} 

= u? - a 1 + -9 2 , (2.27) 
3 

1 1 1 3 

III(u l .j) = -Ui,jUj,kUk,i - -^{Ui,i){Uj,kUk,j) + g(«i,t) 

= ^V r • j~V r • (u(V r • U) - (U • V r )u) U 

- (u(V r • u) - (u • V r )u) • V r u} 



( _ , 20[er +-u ) - <Jijcrj k (T ki - (TijUJiUj (2.28) 
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where we have used the following relations 



' .2 , _2 1/j2 



u i jUij = ef + (T* + -6 A , (2.29) 

^ ».,,».,.. = -u 2 + <? 2 + ^9 2 ■ (2.30) 
Using those invariants, we can describe Q-p as 

Q v = 2(ll(u^)) v - 2 -{l{u^)) 2 v . (2.31) 

In the following discussion, we will use the Lagrangian description. Then, it will be convenient to rewrite 
Eq. ( |2.31 ) by the Lagrangian variables. The Lagrangian domain is just the initial domain T>i. The spatial 



average of a tensor field A is described by the Lagrangian coordinates as 
where the Lagrangian average {A) v is defined by 



d 3 q A (2.33) 



Using this relation, the backreaction term is described by the volume average in the Lagrangian domain X>j 
as 

QV= (^ <JrII(U,j)) ^ _ 5{(J^ (JrIKj)) ^} ' (2 ' 34) 

So far, we have not made any approximation. However, when we evaluate those invariants explicitly, we 
need a further ansatz. Here we adopt a perturbation method. In the conventional perturbation approach, 
we expand the variables around the background Friedmann universe. In fact, Buchert, Kerscher, and Sicka 
evaluated the backreaction term Qx> by this perturbation method [ fl8| |, which we will show in the next 
section. We have also another possibility to divide a perturbed part from a uniform background part as 
follows: Since we average the variables in some domain and take into account its backreaction effect, we 
can adopt the averaged variables as unperturbed parts. This way to extract the perturbed parts seems to 
be more consistent and may provide some difference from the previous approach by Buchert et al. We will 
discuss the detail next. 



III. LAGRANGIAN PERTURBATION EQUATIONS WITH BACKREACTION 

The conventional Lagrangian approximation (for example, ZA 0) was constructed in a homogeneous and 
isotropic universe. However, when a nonlinear structure is formed, it is not so clear what is the background 
universe and how to divide the perturbed part from the unperturbed one. Although the simplest way is the 
conventional approach, as we discussed in previous section, a local expansion rate of the universe in some 
finite domain will be affected by its inhomogeneity. Then the growth rate of density fluctuations may also 
be affected by the inhomogeneity. 

We shall derive the consistent equations for the Lagrangian perturbations with the backreaction term Qv- 
First, we introduce the "comoving" Eulerian coordinates x& and peculiar velocity v& as 

r = a b (t)x b , (3.1) 
u ee r = a b (t)x b + v b . (3.2) 

In the conventional Lagrangian perturbation approach, we set a b = an{t) where an is just a scale factor of 
the background Hubble flow of the whole universe and satisfies the Friedmann equation. However, because 
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of the above reason, the expansion of a domain T> may be described by op. Hence, we have not fixed a scale 
factor a b yet. It will be determined later by adding a further ansatz. 
With the relations 

Vr = —V Xb , (3.3) 

a-b 



the basic equations (p.4|)-(2.7) with P = are rewritten as 



^ + -V Xb .{v b (l + ( 5 b )} = 0, (3.5) 
at a,b 

^f) +— (vb-VxJv 6 + ^v 6 = -— V X6 ($ + la b a b xl) , (3.6) 
at J Xb a b a b a b \ 2 ) 

' 2 2 



A Xb $ - -nGp b aixi = AnGai Pb 5 b - Aa' b , (3.7) 



where the density p is divided into an unperturbed uniform part p b (t), which is defined by the mass conser- 

3. 
}>- 



vation p b a^= constant, and its fluctuation S defined by 



6 b = P —J^. (3.8) 
Pb 

$ is the gravi tatio nal p oten tial defined by g = — Vr^. 

From Eqs. (3.6) and (3.7), we obtain the following equation: 

3^+4tgJ+-V +-(v6-V x Jv 6 +^v 6 l+47rGp6^ = 0. (3.9) 

a b J a b {\ ot J Xb a b a b j 

If we specify a b , this equation provides the equation for inhomogeneous part subtracted a background uniform 
expansion. For example, in the conventional Lagrangian approximation, a b is taken to be a scale factor of 
the FRW universe an, which satisfies the Friedmann equation 

3— + iwGpn - A = , (3.10) 
a H 

where P h is the mean energy density of the universe. On t he ot her hand, if a b is chosen to be ap, which 
evolution is given by the averaged Raychaudhuri's equation (2.21), the perturbation equations are modified 
because of the backreaction effect Qx> as 

— V **> • I (nsr) + — (v© • V^Jvp + I + AkG Pv & v = -Q v , (3.11) 

where the variables with subscript T> are defined those with the scale factor a-p and p-p = {p)x>- 
The rotation of (3.6) yields 

Vx6X {(^) Xb + i (Vb - VxJVb+ S V6 }^°- (3 - 12) 

Now we will describe those equations in terms of the Lagrangian coordinates. Using the Lagrangian time 
derivative along the fluid flow 

d^l + ^ v - v ^ (3 - 13) 
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Eqs. ( 3.11 ) and ( 3.12Q are rewritten as 



- V x, • + -V6 ) + ^G Pb 5 b = -Q b , (3.14) 



a b \ dt a b 

< 3 - i5 > 

where we have introduced Qb, which is zero for the conventional Lagrangian approach without backreaction 
but is chosen to be Qj> when we include the backreaction due to inhomogeneity. Note that there are still 
some space to estimate this backreaction term depending on a choice of the background scale factor a b . 

The Lagrangian coordinates qf,, which follows the background flow, is defined by the initial values of the 
comoving Eulerian coordinates x&, and the Lagrangian perturbations are described as 

x b = q 6 + S b (q 6 ,i) . (3.16) 

where S b denotes the displacement from the uniform distribution assuming that the scale factor is given by 
a b . The continuity equation ( |3.5| ) or equivalently the mass conservation yields 

dM = pd 3 r = pald 3 x b — palJ b d 3 q b — Pid 3 q b = constant, (3-17) 

where J b = det(cfeq,j/9q&j) = det(<5jj + dS b i/dq b j) is the Jacobian of the coordinate transformation x b — > q^. 
With p b al= constant, we have 



or equivalently for density contrast 



p oc p b J b 1 , (3.18) 



J, 1 -I. (3.19) 



Now all physical quantities are found to be written in terms ofS^ and it rem ains only to find solutions 
for St. We obtain the following equations for S b from equations (3.14) and ( |3 .15 ): 

V X6 • (s b + 2^S fc ") - -4irG Pb ( J,- 1 -l)-Q b . (3.20) 



2^S b 

a b 



0. 



(3.21) 



In what follows, the subscript b will be dr opped exce pt for a b and p b , but we may put it when we have 
to distinguish them. In order to write Eqs. (3.20) and ( [3.21 ) in terms of q derivative, we use the relation 



between V x and V q = <9/<9q, which is given by the definition as 



d_ 



dxj d 
dqi dxj 



6 .i> 



dSi 



d 

dx-j 



d 

dxi 



dSj d 
dqi dxj 



(3.22) 



Using this equation iteratively, we 


have 






d 


d 




d 


dxi 


dqi 


dqi 


dx 3 




d 


d^ 


(~ 




dqi 


dqi 


\dqj 




d 


dSi 


d 




dqi 


dqi 


dqj 



dS k d 
dqj dx k 
dSjdS k d 
dqi dq, dx k 



(3.23) 



This gives the differential equations in terms of q formally. In order to write the basic equations explicitly, 
however, we need perturbative approach. We then expand the deviation vector S as S = S*- 1 -* + S( 2 ) + ■ • •. 
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Superscripts (1) and (2) denote first and second order quantities in a perturbative expa nsion with respect 
to magnitude of small displacement e from uniform distribution, respectively. From Eq. ( [2.22 ), the leading 
terms of the backreaction is se cond order of e. Then we have to derive the perturbation equations until at 
least second order. From Eq. ( [3. 15 ), we obtain the first and second order equations as follows: 



V, x (S( 1 )+2^S« 

a b 



V. x I S( 2 >+2^S( 2 > 
a b 



, 



s ijk^>i\j \ Z k \l + Z ~ a k\l 



(3.24) 
(3.25) 



Next we consider equation (|3.14| ). The Jacobian J is expanded as 



J = 1 + V q ■ S« + V q ■ + \ [(V, • s^f s^s% 



+ 



(3.26) 



Thus we obtain the first order equation, 



V.. (s( 1 )+2^S«-4 7 rG /9h S« 



(3.27) 



and the second order one, 



V,. (s^+2^S^-47rGp b S^ 
a b 



(3.28) 



In order to solve the above perturbation equations, we decompose S*- 1 ^ and S^ 2 ) into the longitudinal and 
the transverse parts in the form 



s (2) = v 9 c + C T , v g -c T = o, 



(3.29) 
(3.30) 



where ip and £ are first-order and second-order scalar functions, respectively. As for their physical meanings, 
the first-order longitudinal and transverse parts are related to linear density and vortical perturbations, 
respectively. For the second-order level, however, such a simple interpretation of the perturbation modes 
does not hold any more. 



The first-order perturbation equations (3.27) and (3.24) become 



,, ( ( • I- A\) - \-Cp h i: 

V„ ( «/•' • 

1 a b 



0, 

: . 



(3.31) 
(3.32) 



The second-order perturbation equations (|3.28 ) and ( 3.25| ) are also divided into the longitudinal and the 
transverse parts as 



(3.33) 



, r a b ■ t 



4TTGp b e ijk xpf^^\ k i 



(3.34) 
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In what follows, we will fix ab to be ap. Note that if we set ab = an, the above equations are equivalent to 
those by Buchert et al. Since the displacement vector S denotes a deviation from homogeneous distribution 
in the local domain T> 1 the average density pb — (p)-p, which is written by <5p in what follows, is not equal 
to that in the whole universe. The density fluctuation described by S should be defined as 



P~ Pv 
Pv 



with 







(3.35) 



Here we use the subscript T> to show what is the background mean density. 

To solve the above equations, we carry out further procedure. Since the first order equations (3.31) are 
the same as those of ZA except for a scale factor, we can easily solve them. The perturbation variables are 
separated into time and spatial functions as 



The evolution equation for 6p is 



2 — bj} — 4:irGpT>bT> = 



(3.36) 



(3.37) 



while y>p(q) is determined from the initial data. 

As for the second order equations, we need a little more complicated procedure. In Appendix A, we 
rewrite the backreaction term Qp in context of the Lagrangian approximation and present its explicit form 
up to the second order perturbations (ij),ip T ,(,, and £ T ). Since the transverse modes (tp T and £ T ) may not 
be so important, we shall ignore them in what follows. Then the leading order of the backreaction term Qx> 
becomes very simple as 



(3.38) 



Q B = 2^n(^| ii ))-(i(V'p|i J -; 

where \i = <9 q , (see Eq. ( |A14D . Using ( |3.36[ ), we find 



which depends only on time t. Eq. ( 3.33j ) for £p is now 



A q ( ( v + 2^( v - ^Gp b Cv ) - ^Gp v (^ij^ij ~ (A 9 V) 2 ) - Qb 

&b 



Since the first term of r.h.s. of Eq. ( 3.33| ) 



(3.39) 



(3.40) 



(3.41) 



depends on spatial coordinate q as well as t, we have to divide a solution into two parts: one is inhomogeneous 
term and the other is homogen eous one. First, in order to find a homogeneous term, we commute the time 
and spacial derivatives in Eq. (3.4C), and then integrate it over T>i, finding 



d^ <i-p d 

— (A qv Cv) Vt + ( A 9bCx>>x, 4 - 4:irGp v (A qT> C,v) Vi 

= 2nGp v bp (ip v \ijipT>\ij ~ (A^^x,) 2 )^ - Q v , (3.42) 

which is an ordinary differential equation for (Aq^Cp)^, . By Introduction of new variable (x> by 

A to Cd = A TO Cc ~ (A^Cd)^ , (3-43) 
we can eliminate the backreaction term Qx>. The equation for (p is now separable by setting 
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Cl> = C V (t)xT>(<l) , 



(3.44) 



c v + 2— cp - AnGp v {c v + Ab 2 v ) = , 
clt> 



A 



(3.45) 
(3.46) 



where A is a separa tion c onstant. Finally we obtain a set of basic equations ( 2.21[ ),(3.37), ( |3.42| ), and (3.45) 
with the definition (3.39) for dynamical variables ap, &p, (A^^p}— , and cp. The spatial inhomogeneities 
ip-p (q) and Xx>(l) are determined by initial distributions and Eq. (3.46). 



IV. RELATION BETWEEN TWO LAGRANGIAN DESCRIPTIONS 
: SETTING UP INITIAL DATA 

In order to solve our basic equations for Lagrangian perturbations, we have to set up our initial data. 
Since the initial fluctuation is given by deviation from uniform distribution of the whole universe, which is 
described by the Eulerian coordinates, we have to construct our initial data from those data by an appropriate 
transformation. Furthermore, since our equations are valid only in each domain T>, if we wish to analyze 
some statistical properties of our results, we have to go back to the whole universe. Here we will first discuss 
such a transformation, and then set up our initial data. 

Suppose that we have some inhomogeneous distribution in the Eulerian coordinates r. We then find an 
and dp by integrating a volume over the whole universe and over a domain T>, respectively. Here we set 
these initial values equal to unity. These scale factors fix x# (then q# and S# ), and xp (then qp and Sx>). 
We also find pu and pp by averaging the density p(r, i) over the whole universe and over the domain D, 
respectively. 

The density perturbation in our basic equations is deviation from a density averaged in a domain T>. 
Hence a perturbation in the conventional Lagrangian approximation does not satisfy the condition ( 3.3 5| ) . 



We have to reconstruct initial conditions in our local domain from fluctuations given in the whole universe. 
In order to set up initial conditions, we may need only the relation between the Eulerian coordinates r and 
our Lagrangian coordinate qp. However, we also need the relation between two Lagrangian coordinates, qp 
and q# from the following reason. In an inhomogeneous universe, we do not know in which domain we are 
living. Then, when we analyze our results, we need a statistical analysis in the whole universe to know how 
the observed values are plausible. In particular, one of the most interesting observed values is a peculiar 
velocity, for which we have to introduce a Hubble flow an and the comoving Eulerian coordinate x#. As 
a result, we need a relation between q# and qp as well. Note that although both Lagrangian coordinates 
follow matter fluid, those scales are different. Because we have mass conservation for infinitely small domain 
(qz>, qc + Aqp), or equivalently (q#, q H + Aq ff ), i.e. 

AM = pHa 3 H (Aq H ) 3 = ppa|,(A( ? p) 3 . (4.1) 

Assuming 

av{U) = a H (U) = 1 , (4.2) 

we find 



Aq H 



Ph(U) i1/3 



pv(U) 



(4.3) 



Since the density perturbation in the conventional Lagrangian approximation is given by 



6h = P -^, (4.4) 

PH 
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we find px> = Pb(1 + (£ff)p) from the definition of px>- This gives the relation between g-p and qn as 

We then find the relation between (5// and St> as follows: From 

iii P P PH , x s Ph(M) avtt) 

Pv Ph Pv Pv(U) a H (t) 



with Eq.(|4.5|), we find 



3© 



(4,(1 + 6 B ) 



1 



4 [i + (M*i)>2> 4 j 
We also obtain a peculiar velocity against the Hubble flow as 

vh = r - Hr = cihSh = iflv - Hav) xp + apSp . 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



We will use this definition in analysis of a peculiar velocity. 

mditioi 

5v(U) 



Now we can set up our initial condition. As for the density perturbations, we find from Eq. (4.7) as 

(i + M«0) 



1 + {Sh(U))i 



1 



(4.9) 



From this equation with definitions (3.19), (3.26), (3.29), and (3.36), we find the initial condition for ipx> as 

(6 H (U)) Vz -5 H (U) 



1 + 5 H {U) 



(4.10) 



where we set b-p(ti) — 1. 

As for the scale factor a-p of the do main V, its initial value is set to be unity (Eq. (|4.2D ), but its time 
derivative is non-trivial. In fact, From ( |2.16| ), we find ax> as 



av{U) = a H (ti) 



1+3 (I 



(4.11) 



V. EFFECT OF INHOMOGENEITY 



Here, we study a simple model to show new aspect in our approach. We assume a plane-symmetric 
1-dimensional model. The Lagrangian perturbation is given by 

V^ = S(q)=&(t)(a(gi),0,0). (5.1) 

In the conventional Lagrangian approximation, ZA gives an exact solution in a plane-symmetric case. How- 
ever, it does not take into account a backreaction effect of inhomogeneity on the Hubble expansion. Since 
our approach includes the backreaction effect, we will analyze this simple one-dimensional model and com- 
pare our results with those by ZA. We also look at a difference of the backreaction term estimated by the 
conventional Lagrangian approach (Buchert et al). 

When we consider our backreaction effect, we have a set of equations: 

ij v + + ^GpviPv = -Q-d , (5.2) 

2 (bl(h) v ) 2 

Qv = ~ J —2- (5.3) 

3(1+62 p. 
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Since we analyze the present model numerically, we have to introduce the size of "whole universe", which 
is L. We the introduce the scale of averaging domain I. As for initial conditions, we adopt a power law 
spectrum with the index n = 1: 

Pi{k)(xk. (5.4) 

We also introduce a cutoff at small scale, which wave number is k cu t — 1024/co, where k = 2tt/L. We set 
that the initial time is a = 1. Then the amplitude of fluctuation is chosen so that the first shell-crossing 
occurs at a ~ 1000. The number of grids is N — 2 16 and we use a periodic boundary condition. We take an 
ensemble average over 500 samples, which initial conditions are given by random Gaussian. In our approach, 
since we do not know where we are living, we study its statistical properties. In particular, we will see 
the scale dependence of the averaged variables and the probability distribution of the Hubble parameter, 
deceleration parameter, and pair- wise velocity. 



A. Hubble parameter 



First we analyze the expansion rate of local domain. If we fix the Hubble parameter H Q by local observation 
in a domain T>, the most probable value of H is given by the averaged expansion rate of the domain, which 
is (8) v /3. Several authors so far discussed such a local measurement of the Hubble parameter ]l9|-|2^]. 
In particular, Shi and Turner |23[ estimated a possible value of the Hubble constant measured locally and 
discussed a deviation from the global value, using linear perturbation theory with the CDM model. They 
found that for small samples of objects that only extend to 10,000 km s _1 , the variance can reach 4%, while 
for large samples of ob jects to 40,000 km s _1 , the variance is about 1-2 %. 

We solve Eq. (2.21) for a-p with a backreaction due to inhomogeneity. If we are living in an undcrdense 
region on average, the expansion rate will be faster than the Hubble one for the whole universe. While, if 
we stay in an overdense region, the rate will be slower than the global Hubble one. Fig.0 shows the PDF of 
a local Hubble parameter in our calculation. If our domain is small, the deviation from Hq gets large. For 
example, the dispersion of the Hubble parameter is about 1.2 % for the I = 128-grid domain, while 0.66 % 
for the? = 256-grid domain. The dispersion of our model is consistent with the result by Shi and Turner. 



o.i 
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FIG. 1. The PDF of the Hubble parameter at a = 900 for I = 128-grid domain. The dispersion is 1.2 %. 



B. Density fluctuation 



Next we show the PDF of density fluctuations. In the Eulerian linear approximation, if initial data is 
given by random Gaussian distribution, the PDF of density fluctuations will remain its Gaussian form during 



12 



evolution. On the other hand, in the Lagrangian approximation, there appears a nonlinear effect. In fact, 
Kofman et al shows that the PDF approaches to a log-normal function rather than a Gaussian function in 
the cases of the Lagrangian approximation and N-body simulation p4j . Padmanabhan and Subramanian 
also discussed the PDF with the ZA and found a non-Gaussian distribution |2q| . 

Here we analyze the PDF of density fluctuations using our approximation. The results are shown in Fig.^. 
From comparison with the result of ZA, the void region (i.e. an underdense region; S < 0) is found in higher 
probability in our approximation. Especially, if the size of a domain is smaller, the difference gets larger. 
On the other hand, the probability to find an overdense region (5 > 0) decreases in our approximation. 



-1=128 grids 
ZA 



-1=128 grids 
ZA 




0.1 
0.01 
0.001 
0.0001 
10" ' 



-0.6 -0.4 

s 

(b) 



FIG. 2. The PDF of density fluctuation at a = 900 for I = 128-grid domain and that with the ZA. The 
probability to find an overdense region for I = 128-grid domain is less than that with the ZA((a)). On the other 
hand, the probability for an underdense region increases for I = 128-grid domain((b)). 



The reason is very simple: During evolution, an overdense region shrinks and a nonlinear structure is 
formed as the Zel'dovich's pancake. On the other hand, an underdense region expands. Therefore, although 
the initial volumes of overdense underdense regions are the same, the volume of the latter gets larger than 
that of the former in a nonlinear stage. In addition to this Lagrangian nonlinear effect, we take into account a 
backreaction effect. This effect enhances expansion of an underdense region and contraction of an overdense 
region. As a result, the above difference between ZA and our approximation appears. 



C. Peculiar velocity distribution 



In the conventional Lagrangian approximation, if an initial condition is given by a random Gaussian 
distribution, the PDF of a peculiar velocity also remains its Gaussian form [g4| . Because a peculiar velocity 
vh in the conventional Lagrangian approximation (or the ZA) is given by 

= a H S H , (5.5) 

the spectrum of a peculiar velocity is proportional to that of a density fluctuation as P v (k) oc (b/B) 2 Ps(k). 
However, in our case, a peculiar velocity is given by ( [l.g| ) and the growing factor bx> is different in each 
domain, the PDF could deviate from a Gaussian distribution. However, from our numerical analysis (see 
Fig. ||), the PDF of a peculiar velocity seems to still be a Gaussian. This may be because we need only a 
linear term of the Lagrangian perturbations in the present analysis. 
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FIG. 3. The PDF of peculiar velocity at a = 900 for I = 128 grid domain. The PDF obeys Gaussian form. 

D. Pairwise peculiar velocity distribution 

The PDF of a radial pairwise peculiar velocity is known to show an exponential form from analysis of 
N-body simulation and the ZA [£6]]27]]. For this variable, even if initial data is given by a random Gaussian 
distribution, the PDF approaches an exponential form as the universe evolves ETj ]. The origin of this result 
could be understood by nonlinearity of gravity. In the case of one-dimensional plane-symmetric system in the 
conventional Lagrangian approximation (in fact, the ZA is exact), we do not find any non-Gaussian structure. 
However, even in the ZA, the PDF shows non-Gaussian behavior in the case of three-dimensional case p7j . 
If we take into account our backreaction effect, does non-Gaussian behavior emerge even in one-dimensional 
case ? 

The pairwise peculiar velocity is defined as follows: 

VAs(i) = v\b(£) - v A (t) 

= vy(i) + vj_, (5.6) 

where V|| and vj_ represent components parallel and perpendicular to xab = — x^, respectively. In a 
plane-symmetric case, V|| only appears. Hereafter we write this by v. If matter distribution is clustering, v 
is expected to be negative. 

Giving initial data by random Gaussian, the PDF of a pairwise peculiar velocity is Gaussian at initial 
time. During evolution, the PDF will deviate from Gaussian. In fact, from Fig.[|, which shows the PDF of a 
peculiar velocity in nonlinear regime, we find that it is not Gaussian and approaches an exponential form in 
a small velocity region. As a reference, we show that the PDF for the ZA, which shows a Gaussian form. The 
reason may be understood as follows: In a plane-symmetric model, gravitational potential is proportional 
to a distance between two sheets (ip ~ \r\). Then, even if two sheets approaches very closely, a gravitational 
force does not become strong but keep constant. When we take into account a backreaction effect, however, 
a gravitational force will be strong in a clustering region, because the expansion rate of a local domain is 
slow down. The strengthening of a gravitational force in a cluster region may make a deviation of the PDF 
of pairwise peculiar velocity from its Gaussian form. Note that in the 3D system, which shows non-Gaussian 
PDF, a gravitational force increases when two particles approach. 
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FIG. 4. The PDF of pairwise peculiar velocity at a = 900 for I — 128-grid domain (a) and E-dS model (b). The 
PDF in (a) does not obey Gaussian form, and approaches to log-normal form at small velocity region. On the other 
hand, the PDF in (b) obeys Gaussian form. 



E. Deceleration parameter q 



Another interesting observable variable is a deceleration parameter. The recent observation of type la 
supernova may suggest an acceleration of the Universe (29) . Although this result may naively suggest an 
existence of dark energy such as a cosmological constant A, we could find some effective model without dark 
energy which explain the observation. Then we shall estimate a deceleration parameter averaged in a local 
domain here. 

We define local deceleration parameter qx> as 



qv 



avav 



"7 



(5.7) 



which can be evaluated by Eq. ( 2.21 ) and -ffp- Buchert et al showed the evolution of deceleration 
parameter for the ZA. They picked up overdense and underdense regions of three-cr fluctuations and found 
that qx> for an underdense region could be a present day value, which is smaller by more than 200% than 
that of the background E-dS Universe, although such a region is still decelerating. 

Although our approach includes a backreaction consistently, our analysis shows that a deviation of qx> 
does not get so large. The difference of q-p from the ZA is very little even just before the shell crossing. We 
show the time evolution of q-p for a plane-symmetric 1-dimensional model in Fig. 0. Even if a domain is 
extremely underdense, the domain is decelerating. This may be because our approach is still perturbative. 
We will discuss it further in the next section. 
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FIG. 5. The evolution of a deceleration parameter qx> for I — 64-grid domain. The overdense domain shows 
increase of qv (a solid line for maximum value), while the underdense domain shows decrease of qx> (a dashed line 
for minimum value). 



VI. CONCLUSION AND REMARKS 



We propose new Lagrangian perturbation theory with a backreaction effect by inhomogeneity of density 
perturbations and present a set of basic equations. The inhomogeneity affects the expansion rate in a 
local domain and its own growing rate. In a one-dimensional plane-symmetric model, we have numerically 
analyzed our approach, and calculated the growing rate density perturbations and the PDF of several 
observed variables. We set our initial conditions as random Gaussian distribution. From our analysis, we 
show that the expansion rate of an overdense region is faster than that of the whole universe as expected. We 
also show that the local Hubble parameter may deviate from the global one by about 1.2 % for a I = 128-grid 
domain. It may be too small to distinguish its effect in the present observations ]2q] , but it could become 
important in fully non-linear stage, which we cannot describe in the present approach. 

The PDF of density is slightly different from that of the ZA. In our model, an underdense region expands 
faster than that in the E-dS model and then its volume gets larger. Hence, a probability for a negative 5 
region increases as seen in the PDF. As for a peculiar velocity, even if we take into account a backreaction, 
its PDF is still Gaussian. 

The PDF of pairwise peculiar velocity, however, shows an effective difference from the conventional La- 
grangian approach. In one-dimensional plane symmetric case, the PDF in the conventional Lagrangian 
approximation (the ZA) is Gaussian, but ours is not but approaches an exponential form in a small relative- 
velocity region, which agree with the N-body simulation and the 3D Lagrangian approximation f26|j27|| . 

Finally, we mention about recent observation about cosmological parameters. According to the observation 
of type la supernova, the expansion of the Universe seems to accelerate 1 29 . Combining observation of the 



cosmic microwave background radiation (CMBR), the result suggests existence of dark energy such as a 
cosmological constant A [ p0| . However, this produces another difficulty, that is the so-called cosmological 
constant problem. To avoid such a difficulty, if we could explain the observation without cosmological 
constant, it would be more natural. Recently, Tomita discussed such possibility assuming we are in a large 
local void |3l|-|33T]. Globally the Universe is flat (EdS universe), but we are sitting near the center of a 
local void, which existence is observationally confirmed. Then he calculated the luminosity distance, finding 
that the observation can be explain by such a model. We then wonder whether we could have the similar 
explanation if we are living in an effective void, which is a domain with an averaged energy density below 
the critical value. Since we have to treat a strongly nonlinear structure to explain the observation, this 
is beyond our present approach. However there is some indication. If we analyze the time dependence of 
the backreaction term Qv, which evolves as a -1 in a linear perturbation level. This time dependence is 
the same as a perfect fluid with the equation of state P — — |p, which could be dark energy. If we could 
explain the above observation without introduction of any strange matter but just by inhomogeneity of 
density distribution in the Universe, we will have a natural understanding of the Universe. This is under 
investigation. 
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In addition, there are two further approaches which are related to the present work. Takada and Futamase 
jjij proposed that they divided Lagrangian perturbation to large-scale and small-scale perturbations, then 
discussed interaction between those scales. Taruya and Soda |35| discussed dynamics of averaged variables 
in the case of a spherical infall model taking into account a backreaction effect. Since those are interesting 
approaches, it may be useful to use their approaches to discuss the present subject and compare those results 
with ours in future. 
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APPENDIX A: BACKREACTION IN LAGRANGIAN APPROXIMATION 



Using the Lagrangian approximation, we estimate a backreaction term Qx> ■ Using the Lagrangian pertur- 
bation S, we shall rewrite the backreaction term. For convenience, first we define the functional determinant 
of three functions A(q), B(q), C(q): 



1(A, B, C) = df ' -1 9- _ SijkA^C^ . 
The invariants of the velocity gradient Ui j are given in terms of functional determinants as 



1 
1 



(Al) 

(A2) 
(A3) 
(A4) 



Introducing The r.h.s. of the Lagrangian coordinates qi and perturbations Si by r» = a(t)(qi + Si), the r.h.s. 
of these equations are rewritten as 



J(h,rj,r k ) = a 2 a [J(qi, qj,qk) + 3J(Si,qj,q k ) + 3J(Si, Sj, S k ) + J(Si,Sj,Sk)] 



-a 



2 

-2a 2 a 
-a 3 



J(Si, qj ,q k ) + 2J(S t , Sj,q k ) + J (Si, Sj,S k )\ , 
J(h,fj,r k ) = ad 2 [J(qi,q 3 ,q k ) + 3J{S l ,q j ,q k ) + 3J(Si, Sj,q k ) + J (Si, Sj, S k )] 
J(S h q,,q k ) + 2J(S h Sj,q k ) + J(S U Sj,S k ) 

J(Si,Sj,q k ) + J(Si,Sj,S k ) . 
Introducing a simplified description follows: 

h = i(%) , n s = ii(%) , = m(%) , 
h = i(%) , lis = n(%) , nii. = m(%) , 

we find that most terms in (|A2])-(|a"4|) with (|A5|), (|A6|) are given by the following quantities; 



Sij k J(qi,qj,q k ) = 6 , e ljk J(S l ,q j ,q k ) = 2I S , 
£ijkJ(Si,Sj,q k ) = 211 s , £ij k J(Si, Sj, Sk) = 6III S , 

£ij k J(Si,qj,q k ) = 21 s , £ij k J(Si, Sj,q k ) = 211 s . 



(A5) 



(A6) 



(A7) 



(A8) 



J (Si, Sj,q k ) and J (Si, Sj, S k ) are written as 



£i 3 kJ(S l ,S j ,q k ) = —II, , 
d 

£ijkJ(Si,Sj,S k ) = 2— III S , 

but the last term J (Si, Sj, S k ) cannot be written by any of I s , II S , III S and those time derivatives. 
We then finally obtain 



(A9) 
(A10) 



1 



J r l( Uiij ) = - i 6a 2 a (1 + I s + II S + III S ) + 2a A I, + 2— II S + —III 







d 



dt 



dt 



(All) 



J r II(u id ) = - [6aa 2 (1 + I s + II S + 6III S ) 
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+ 4^(l + I s + |lI s + |lIL 
+2a 3 [lh+e ijk J (Si, Sj.S*)}] . 



(A12) 



Using these invariants, we find a backreaction term in terms of the Lagrangian perturbations as 

2 1 

Q-d = - 



3 (1 + (I s ) Vi + (U s ) Vt + (llh) Vt ) 2 
x [3(1 + (I s ) Vi + (II s ) Vi + (III.)^) 

(Hj}^ + (^ £ ijkJ(Si, Sj, Sk)^ ) 



(A13) 



If we take into account the perturbations up to the second order, we find 



Qv = 2 (II&) Vi - - (lv,^)p. , 



(A14) 



where we use the variables in Eqs. ( |3.29 ) and (3.3C). Note that the second order perturbation £ do not 
appear in Eq. ( A14 ). 

Buchert et al estimated the backreaction term Qx> assuming ZA, i.e. a?, = a# and the Lagrangian 
perturbation S is separable as 



x = q + 6(i)V g ^(q) 



(A15) 



and found the similar form to (A13). flli 
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